c
c******************************************************************
        function pm(xmu)
c****************************************************************
c
       implicit real*8(a-h,o-z)
       pm    =0.5233-0.0815*xmu+.0135*xmu*xmu-.001*xmu*xmu*xmu
       return
       end
